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Abstract. - We examine the global organization of growing networks in which a new vertex 
is attached to already existing ones with a probability depending on their age. We find that 
the network is infinite- or finite-dimensional depending on whether the attachment probability 
decays slower or faster than (age) -1 . The network becomes one-dimensional when the attachment 
probability decays faster than (age) -2 . We describe structural characteristics of these phases and 
transitions between them. 



A relation between small worlds (compact, infinite- 
dimensional objects) and large worlds (finite-dimensional 
objects) is a key issue in understanding of networks. The 
surge of interest in networks has been triggered by the pa- 
per of Watts and Strogatz [1] who found a transition from 
lattices to networks with a small-world organization. The 
Watts-Strogatz construction describes a smooth crossover 
between the two limiting situations. A sharp transition 
between finite- and infinite-dimensional networks (it co- 
incides with a so-called searchability point) has been ob- 
served in studies of decentralized search algorithms [2-4]. 
These findings have been made in "static" networks. Re- 
sults of numerical simulations [5-8] suggest a sharp tran- 
sition also in specific growing network models. The exis- 
tence of this special point was confirmed analytically in 
work [9]. Yet essentially nothing is known about the na- 
ture of these transitions in static and especially in growing 
random networks. Here we provide an analytical descrip- 
tion of a transition from an infinite-dimensional network to 
a finite-dimensional one in random recursive graphs with 
aging. 

A random recursive graph is one of the two principal 
models of a random graph (the second one is the Erdos- 
Renyi model). A random recursive tree (a tree is a con- 
nected graph without loops) is built according to the fol- 
lowing procedure — a new vertex attaches equi-probably 
to an already existing vertex. Random recursive trees are 
well understood, see e.g. Refs. [10-16]. Here we generalize 
this model by letting the choice of the target vertex to be 



proportional to a certain function, w{a), of the age. (The 
vertices in the network are labelled as i = 0, 1, 2, . . . , N, 
where i = is the oldest vertex. So the age of vertex i 
is a = N — i.) Thus the probability that vertex N + 1 
attaches to a vertex of age a is Pat (a) = w(a)/W(N), 
where W(N) = J2i <a<N w(a). To be more specific, we 
shall consider attachment kernels with an algebraic large 
age asymptotics, 



w(a) ~ Ao 



when d> 1. 



(1) 



For this w(a), the normalization factor W(N) converges 
when fi > 1 and diverges otherwise. Degree distributions 
in scale- free models of this kind were analysed in Ref. [17]. 
As exponent fj, decreases from +oo to — ao, the resulting 
network transforms from a chain to an ultimately com- 
pact star graph. Here we outline how the network varies 
with exponent p. Details of derivations will be presented 
elsewhere [18]. 

Results. — The simplest geometric characteristic of 
the network is the mean distance £(N) of the Nth vertex 
from the root. This quantity behaves similarly to the mean 
intervertex distance and to the diameter of the graph. We 
find that the mean depth of the last vertex scales as 
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where W\ = £ a>1 aw ( a ) is the first moment of the at- 
tachment rate [W± is finite when /j, > 2]. The amplitudes 
Ci(/x) and C2(//) in Eq. (2) are 



1 



7+^(2-M)' 



c 2 (a*) 



/i— 1 sin(7r/i) 



(3) 

where 7 = 0.5772 ... is Euler's constant and ip(x) is the 
digamma function, ip(x) = T'(x)/T(x). So, for \i < 2, the 
asymptotics of £(N) is determined solely by the asymp- 
totics of w(a). Remarkably, if fi < 1, the result depends 
only on /i. The mean depth of vertices in the network, (£), 
varies with iV similarly to £(N), Eq. (2), only correspond- 
ing coefficients differ when [i > 1. 

Equation (2) shows that there are three different 
regimes: (i) when fi < 1, the network is infinite- 
dimensional; (ii) when 1 < /1 < 2, the network is rf- 
dimensional, c? = 1/(^—1) (d is the Hausdorff dimension); 
(iii) when fj, > 2, the network is a one-dimensional ob- 
ject. The drastic changes take place at /j, = 1 and /i = 2, 
when the zeroth moment Wo = W = £a>i w ( a ) an d Wi 
become divergent, respectively. 

The mean depth £(iV) is a basic measure of the size of 
the network. A more detailed characteristic is the depth 
distribution Q/{N), that is the mean number of vertices 
at distance £ from the root. We find a set of qualitatively 
different behaviors, and an additional special value /1 = 3 
where the second moment W2 — J2 a >i a2w ( a ) begins to 
diverge, demarcates different regimes. The major transi- 
tions [19] in the behavior of the depth distribution still 
occur at yu = 1 and fi = 2: 

(i) when fi > 2, this distribution is essentially a step 
function with a smeared front evolving as £(N); 

(ii) when fi < 1, the distribution has a well-defined peak 
at - £(N). 



The width of the front or peak scales as 



width oc < 
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We also probed the degree distribution rife. The de- 
pendence of this quantity on the exponent fi is less pro- 
nounced, namely the degree distribution exhibits an expo- 
nential decay when fi < 1 and a factorial decay (oc 1/fc!) 
when fj, > 1. The degree distribution has a surprising be- 
havior at the point /j, = 1 where the transition from an 
infinite- to a finite-dimensional geometry occurs. At this 
point 

n k = 71 TTf 6 " 1 - ( 5 ) 



(fc-1)! 



Intriguingly, the same degree distribution characterizes an 
equilibrium statistical ensemble of random trees [20] . 

Mean depth. — In recursive trees, the addition of 
new vertices does not change the distances from the root 
of already existing vertices. This feature makes the prob- 
lem analytically tractable. Consider first the general case 
when the attachment probability is proportional to an ar- 
bitrary function w(i,k) of the birth time k of a new vertex 
and of the birth time i of an already existing vertex. Then 
the probability that vertex j is at distance £ from 

the root satisfies 



Mo) 



1 



(6) 



n o(j) = $o,j, where the root vertex is j = 0. Consequently 
for the average values of the integer powers m of depths 
£i of individual vertices, we have 



((ii in = 



1 



-5>(UXC>- (7) 



Below we will discuss only m = 1 and use the notation 
£{i) = (£i) for the mean depth of vertex i. In our networks 
w(i,j) — w(j — i) and Eq. (7) gives 



£{N + 1) = 1 + 



W(N)^ 



^Tw(a)£(N + l-a). (8) 



We do not solve this equation directly. Instead, we first 
guess a form of the solution and then verify its consistency. 
We start with a long- memory regime fi < 1. We know that 
for classical random recursive trees (/1 = 0), the mean 
depth grows logarithmically with network size. Suppose 
that the same is valid when \i < 1, that is 



£{N)=Ci InJV + £>! + .. 



(9) 



where C\ and D\ are some functions of fi. To verify (9) 
we first re-write recurrence (8) as 



1 = 



W(N) ^ 



N 

E 



w{a) [£(N + 1) - £(N + 1 - a)]. (10) 



We then plug (9) into (10), make the substitutions w(a) — > 
Aa~^ and W(N) -» AN 1 ^^ /(l- fi), and replace the sum- 
mation by integration E^Li — * ^ Jo 1 21 = a /^- After 
these transformations, Eq. (10) reduces to an equation 
that determines the amplitude 



Ci(/x) 



-(1-m) 



ln(l - x) 



(11) 



Computing the integral (11) yields the resulting expres- 
sion (3) for Ci(/x). The correction Di(/j,) in Eq. (9) de- 
pends on the initial condition and therefore it cannot be 
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computed in the realm of the continuous approach. Note 
the asymptotics: d(u -> 1 - 0) ^ (6/tt 2 ) (1 - u)" 1 and 

d^-ooJe^lnM)- 1 . 

In the marginal case (1 = 1, W(N) ^ A In TV. The 

divergence of the amplitude Ci()i) as p, — > 1 — implies 

that in this case £{N) grows faster than In N. An algebraic 

growth is inconsistent with Eq. (10) where W(N) is set 

to In AT, while trying an accelerated logarithmic growth 

£(N) = B(\nN) x we found consistency when 

When /i > 1, recurrence (8) simplifies to 

N 

l = 5>(a)[*(JV)-*(iV-a)] (13) 



in the large N limit. Expanding £(A — a) into a Taylor 
series and keeping only the first term of the expansion we 
recast (13) into a differential equation 



1 



o = l 



(14) 



If the sum on the right-hand side converges (/i > 2), 
we have d£/dN = l/Wi which leads to the announced 
expression of Eq. (2) in the regime u > 2. The sum 
J2i< a <N aw ( a ) diverges as A In N when fi — 2. Solving 
dl/dN = (AlnN)^ 1 yields £(N) = N/(A\nN). 

When 1 < p < 2, the sum on the right-hand side 
of Eq. (14) diverges as N 2 ~^, so the asymptotics of 
the solution must be £(N) oc N^ 1 . Inserting £(n) = 
A~ 1 C2{p) n^ 1 into (13) and replacing the sum by an in- 
tegral we obtain 



1 



c 2 (m) / dxx-» [i - (i - x y- 1 } , 

Jo 



which yields the resulting expression for C2(p) in Eq. (3). 
Note the asymptotics: C2(p — * 2 — 0) = 1/(2 — p) and 
C 2 (p -» 1 + 0) = (tt 2 /6)(u - 1). Compare the latter with 
Cilfj, -» 1 - 0) ^ (6/tt 2 ) (1 - /i)- 1 which we obtained 
above. 



Depth distribution. 

Eq. (6) takes the form: 



In the graph with aging, 
w(n + 1 — i) 



»=i 



W{n) 



■nt(i). 



(15) 



The average number of vertices at the £ th layer is equal to 



AT 

Qt(N)= ]T 

n=£+l 

Equation (15) may be rewritten as 

n-l 



(16) 



w(j) 



W(n- 1) 



(«-.?)- 7T£+l(«)]- (17) 



In the "short memory" regime (u > 1), W(n) — > VK(oo). 
Since u»(j) quickly decays as j increases, the major contri- 
bution to the sum on the right-hand side of Eq. (17) oc- 
curs when j is small. This prompts us to use the expansion 
w t (n-j) = Mn)-jdMn)/dn+{l/2)j 2 d 2 7: e (n)/dn 2 + - ■ ■ 
and allows us to transform (17) into 



d_ 

d£ 



+ 



d \ If d 2 d 2 \ 

Wl -y t{ n) = -{w 9Wa -^y t (n) (18) 

The above derivation is valid when the moments W\ and 
W2 = J2 a >i a2w ( a ) are finite (that is, when p > 3) and £ 
is large. 

To understand the form of the solution, it is sufficient to 
ignore the right-hand side of this equation, which leads to 
tt^ti) = f(£-n/W{), where / dxf(x) = 1. With Eq. (16), 
this results in Q e (N) = Wi[l-8(£- N/W{j\, where 6{x) is 
the Hcaviside step function. Taking into account the right- 
hand side of Eq. (3) shows that this propagating front is, 
actually, smearing. In the front region we change variables 
(n, £) — ► (t = n,x = £—n/Wi), Tre(n) = p(t, x) and arrive 
to a diffusion equation, 



dp 
di 



D 



d 2 p 



dx 2 ' 

Solving this equation leads to 
Wi 

Qe(N) = -± crfc 



D 



_ W 2 - W 2 

~ zw? ■ 

N - Wi£ \ 
W 1 V4JJNJ 



(19) 



(20) 



where erfc(x) = (2/y / 7r) dte~ l is the complementary 
error function. The resulting width of the front is propor- 
tional to y/ADN. This expression gives the front width 
growing as y/~N when p > 3 and also allows us to obtain 
estimates listed in Eq. (4) in the range 1 < p, < 3. In the 
latter region we simply take into account the dependence 
of the diffusion coefficient D, Eq. (19), on p. 

In the long memory regime (p < 1), vertices of all ages 
continue to evolve and the appropriate variables are r = 
Inn, £ = I — Cilnn, iti{n) = p(r, £), with C\ = C\{p) 
given by Eq. (3). In these variables, we obtain a diffusion 
equation, 

dp^DCh^p 

0r~ 2 de' 1 ' 

where D = C\ (1 - p) /„* dx X -»[\n{l - x)} 2 - 1. The so- 
lution of Eq. (21) is a Gaussian which, after substitution 
into Eq. (16), gives 



Qt(N) 



N 



exp 



(£-(£)) 2 



2D{£) 



(22) 



where {£) = d In A. This is a peak of width - (£) 1 ^ 2 and 
not a step, in contrast to the short memory regime. 

Equation (22) remains essentially valid in the marginal 
case — Here (£) — 3 • 7r~ 2 (lniV) 2 , and the divergence 
D oc (1 — p)^ 1 implies that one must replace D by In A. 
Since the latter is proportional to (t) 1 ! 2 we conclude that 
the width of the peak in the marginal case is enhanced 
it scales as (£) 3/4 . 
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Degree distribution. — The degree distribution can 

be written as N/. = J^^Li c fc( a 7 N) = Nn^, where c/j(a, N) 
is the average number of vertices of degree k and age 
a. In the short memory regime, the joint distribution is 
(asymptotically) independent on the network size N , that 
is, Ck(a,N) — > Cfe(o). This simplifies calculations. On the 
other hand, Cfe(o) changes only when the age is small, and 
hence we cannot use a continuum approximation. The 
equation for Cfc(a), 

c fc (a + 1) = c fe (a) + tu(a)[cfc_i(a) - c fe (a)], (23) 

can be solved and it leads [18] to the following universal 
asymptotics of n/, : 

n>k+i fx TT as - * oo. (24) 
fc! 

This asymptotics is valid for the general class of attach- 
ment kernel (1), where fi > 1. 

In the evolving regime (/i < 1), the density Ck(a, N) 
continues to evolve as the network size grows. In this case 
we employ a continuum approximation. The governing 
equation for Ck(a, N) reads 

f d 9 \ w(a) , . 

U + ^vJ Cfe = w^v) (cfe - 1_Cfe) - (25) 

When fj, < 1, the analysis [18] of Eq. (25) gives an expo- 
nential asymptotics of the degree distribution: 

e -7-V>(i-M) fi _ u \ k 
n k = — — £ when k > 1. (26) 

In the marginal case (/i = 1), an inverse factorial form of 
the degree distribution, Eq. (5), is valid for all k. 

Discussion and summary. — Several network fea- 
tures are puzzling and deserve further attention. 

(i) The degree distribution of our growing random graph 
at (i = 1 is the same as for the static ensemble of random 
trees [20] . We have no explanation for this surprising coin- 
cidence. The two networks differ dramatically: our graph 
is a small world at fi = 1, while according to Ref. [21] a tree 
which is randomly chosen from an ensemble of all trees is 
finite-dimensional with a high probability (the Hausdorff 
dimension of such trees is d = 2). Interestingly, at the sec- 
ond transition point \x = 2, the degree distribution does 
not exhibit any special behavior. 

(ii) The asymptotic (£)n ~ (In TV) 2 arises at a search- 
ability point of specific networks discussed in Refs. [2, 3] 
in the context of decentralized search algorithms. These 
networks are (i-dimensional lattices with random shortcuts 
which connect vertex pairs with probability decaying as a 
power-law of an Euclidean separation. The greedy algo- 
rithms work rapidly at some value of the decay exponent 

— the searchability point. At that point, the run time of 
these algorithms is of the order of (£}n- The coincidence 
of our result with that of Refs. [2,3] suggests that the tree 



approximation is valid at the searchability point of (some) 
networks with loops. 

We have studied simplest recursive graphs which have 
rapidly decreasing degree distributions. Yet similar phe- 
nomena should be observed in growing networks with more 
complex (e.g., skewed) degree distributions. Equations (6) 
and (7) allow one to analyse recursive trees with arbitrary 
kernels w(i,j) generating various degree distributions. 

In summary, we have demonstrated a rich set of regimes 
in a global organization of growing networks with aging. 
Our results have been derived in the context of a specific 
model, yet the major conclusions should be applicable to 
a larger class of networks. Our analysis has exploited the 
absence of loops. Accounting for loops is a challenging 
direction for future work. 

* * * 

This work was partially supported by projects 
POCTI: FAT/46241, MAT/46176, FIS/61665 and BIA- 
BCM/62662, and DYSONET. We thank A. V. Goltsev, 
M. Hase, M. E. J. Newman, A. N. Samukhin, and B. Wa- 
claw for useful discussions. 

REFERENCES 

[I] Watts D. J. and Strogatz S. H., Nature, 393 (1998) 
409. 

[2] Kleinberg J. M., Nature, 406 (2000) 845. 

[3] Kleinberg J., Proceedings of the International Congress 

of Mathematicians (ICM), Madrid, August 22-30, 2006, 

edited by M. Sanz-Sole, J. Soria, J. L. Varona and J. 

Verdera (EMS Publishing House, Zurich) 2007, p. 1019. 
[4] Moukarzel C. F. and Argollo de Menezes M., Phys. 

Rev. E, 65 (2002) 056709. 
[5] Holme P., Physica A, 377 (2007) 315. 
[6] Hinczewski M. and Berker A. N., Phys. Rev. E, 73 

(2006) 066126. 

[7] Tian L., Zhu C.-P., Shi D.-N., Gu Z.-M. and Zhou T., 

Phys. Rev. E, 74 (2006) 046103. 
[8] Zhang Z., Zhou S., Wang Z. and Shen Z., J. Phys. A, 40 

(2007) 11863; Zhang Z., Zhou S., Shen Z. and Guan J., 
Physica A, 385 (2007) 765. 

[9] Lambiotte R., J. Stat. Mech., (2007) P02020. 
[10] Smythe R. T. and Mahmoud H., Theor. Probab. Math. 
Statist, 51 (1995) 1. 

[II] Moon J. W., London Mathematics Society Lecture Notes 
Series, Vol. 13 (Cambridge University Press, Cambridge) 
1974, p. 125. 

[12] Szymanski J., Theor. Comp. ScL, 74 (1990) 355. 
[13] Mahmoud H., Probab. Eng. Info. ScL, 5 (1991) 5. 
[14] Dobrow R. P., J. Appl. Prob., 33 (1996) 749. 
[15] Drmota M. and Gittenberger B., Random Struct. 

Alg., 10 (1997) 421; Drmota M. and Hwang H.-K., Adv. 

Appl. Probab., 37 (2005) 321. 
[16] Krapivsky P. L. and Redner S., Phys. Rev. E, 63 

(2001) 066123; Phys. Rev. Lett, 89 (2002) 258703. 
[17] Dorogovtsev S. N. and Mendes J. F. F., Phys. Rev. 

E, 62 (2000) 1842. 



p-4 



Transition from small to large world in growing networks 



[18] Dorogovtsev S. N., Krapivsky P. L. and Mendes 
J. F. F., in preparation. 

[19] These transitions are different from critical phenomena; 
for a discussion of this subtlety in the context of networks, 
see Dorogovtsev S. N., Goltsev A. V. and Mendes 
J. F. F., arXiv.0705.0010 [cond-mat], (2007). 

[20] Bialas P., Burda Z., Jurkiewicz J. and Krzy- 
wicki A., Phys. Rev. E, 67 (2003) 066106. 

[21] Renyi A. and Szekeres G., J. Austral. Math. Soc, 7 
(1967) 497; G. Szekeres, in: Combinatorial Mathemat- 
ics X, 1982, Lecture Notes in Mathematics, Vol. 1036 
(Springer- Verlag, Berlin), 1983, p. 392. 



p-5 



